First order phase transition in a 1+1-dimensional nonequilibrium wetting process 
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A model for nonequilibrium wetting in 1+1 dimensions is introduced. It comprises adsorption and 
desorption processes with a dynamics which generically does not obey detailed balance. Depending 
on the rates of the dynamical processes the wetting transition is either of first or second order. It is 
found that the wet (unbound) and the non-wet (pinned) states coexist and are both thermodynam- 
ically stable in a domain of the dynamical parameters which define the model. This is in contrast 
with equihbrium transitions where coexistence of thermodynamically stable states takes place only 
on the transition hue. 
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Wetting phenomena are observed in a large variety of 
physical systems in which a planar substrate is exposed 
to a gas phase. The interactions between the substrate 
and the molecules of the gas phase lead to the formation 
of a liquid film close to the surface ^ . By changing the 
temperature or the partial pressure of the gas, such sys- 
tems may exhibit a wetting transition from a phase where 
the thickness of the film stays finite to a phase where the 
film grows and eventually reaches a macroscopic size. 

Much is known about equilibrium wetting transitions 
both theoretically and experimentally. Theoretically, 
they have been modeled by considering the binding pro- 
cess of a fluctuating interface to a substrate. By using 
transfer matrix formulation, it has been found that in 
1+1 dimensions, when the interaction between the in- 
terface and the substrate is short range, the depinning 
transition is continuous |0]. However, when algebraically 
decaying long range interactions are taken, the transition 
could become first order 0. On the other hand, in 2+1 
dimensions the wetting transition is generally expected 
to be first order. However, disorder effects, such as those 
existing in porous media, could make the transition con- 
tinuous W. 

When passing to non-equilibrium phenomena, most of 
the efforts have been made in the characterization of the 
wet phase, namely the description of the growth process 
far away from the substrate. Various models have been 
introduced to account for observed scaling behaviour in 
homogeneous as well as disordered environments ||5|]. 

Concerning the transition itself, a simple 1+1 dimen- 
sional model has recently been introduced where a con- 
tinuous transition is observed in the absence of any bind- 
ing force between the interface and the substrate [g[. It 
would be of interest to study the nature of the transition 
as the interaction between the interface and the substrate 
is varied. 

Motivated by the knowledge of the behavior of equi- 
librium systems, we generalize the model introduced in 



Ref. H by adding an attractive binding force between 
substrate and surface layer. We find that even short 
range attractive interaction is capable of making the tran- 
sition first order. In fact, for a suitable choice of the pa- 
rameter values the model satisfies detailed balance and 
thus we can test our numerics against theoretical predic- 
tions in the context of equilibrium wetting. Additionally 
we can explore the wider parameter space to clarify the 
role of the binding force under nonequilibrium conditions. 
The model of Ref. [g[ is characterized by an adsorption 
rate q and a desorption rate p and exhibits a continuous 
wetting transition at a certain threshold qdp)- The tran- 
sition is related to an unpinning process of an interface 
from a substrate, which may be described in general by a 
Kardar-Parisi-Zhang (KPZ) equation in a hard-core po- 
tential 01 . The additional short-range force is introduced 
by modifying the growth rate go at zero height. Thus, for 
Qo < q (90 > q), there is an attractive (repulsive) inter- 
action between the substrate and the bottom layer. We 
find that a sufficiently strong attractive interaction mod- 
ifies the nature of the unbinding transition, making it 
first order. Moreover, we observe that pinned and mov- 
ing states coexist as thermodynamically stable states in 
a whole region of the parameter space rather than on a 
line, as in the case of equilibrium transitions. This kind 
of behavior has been observed in the past in other non- 
equilibrium models ||^,|9|] and should be a generic feature 
of wetting transitions under nonequilibrium conditions. 

A. Definition of the model: The model is defined 
in terms of growth of a one-dimensional interface on a 
lattice of N sites with associated height variables hi = 
0, 1, . . . , 00 and periodic boundary conditions. We con- 
sider a restricted solid-on-solid (RSOS) growth process, 
where the height differences between neighboring sites 
can take only values 0, ±1. In addition, a hard-core wall 
at zero height is introduced. The model depends on three 
parameters q,qa, and p. It evolves by random sequential 



updates, i.e., in each update attempt a site i is randomly 
selected and one of the following processes is carried out: 

- adsorption of an adatom with probability qq At at the 
bottom layer hi = and probability q At at higher 
layers hi > 0: 



1 



(1) 



desorption of an adatom from the edge of a terrace 
with probability 1 At: 



min(/ij_i,/ii,/ii+i) . 



(2) 



desorption of an adatom from the interior of a terrace 
with probability p At: 



1 if /1.-1 



ij+i 



>0. 



(3) 



A process is carried out only if the resulting interface 
height hi is non-negative and does not violate the RSOS 
constraint \hi — hi±i\ < 1. The time increment per sweep 
(N attempted updates) is At < 1/ max(l, qo,q + p). 

The phase diagram for the case qo — q has been stud- 
ied in Iqjllfl, where a continuous wetting transition was 
found. Clearly, the moving state is not affected by go 
and thus the transition line above which it is stable re- 
mains unchanged. However, the stability of the pinned 
state strongly depends on qo, modifying the phase dia- 
gram and the nature of the wetting transition. In order 
to gain some insight into the mechanism leading to first- 
order transition, we first consider the p — 1 case. Here 
detailed balance is obeyed 0, wherefore the transition 
can be described in the framework of equilibrium statis- 
tical mechanics flQ]. We then consider the case p ^ 1 
numerically. 

B. The case p = I: For p = I and q < 1 the dynamic 
rules satisfy detailed balance and the probability of find- 
ing the interface in a configuration {hi, . . . , h^} can be 
expressed in terms of a potential V{h) by 



P{hi,.. .,/lAr) = ^^^cxp 



N 



-Y.^{h,) 



(4) 



where the partition sum Z^ = ^f^ f^ e~ ^i 
runs over all interface configurations obeying the RSOS 
constraint. The potential is given by 



cxD if /i < 

V{h) = { -Hq/qa) if /i = 

-h \n{q) ii h>Q 



(5) 



As shown in the inset of Fig. H, the attractive interaction 
between substrate and bottom layer is incorporated as a 
potential well at zero height. For q < 1 the slope of the 
potential is positive so that the interface is always pinned 
to the wall. For q > 1, where the slope is negative, 
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FIG. 1. Phase diagram for p = 1. The discontinuous 
(continuous) part of the transition line is represented by a 
solid (dashed) line. The inset illustrates the potential V(h) 
inEq. (5). 

the interface can 'tunnel' through the potential barrier 
and eventually detaches from the substrate. It should be 
noted that in this case, the equilibrium distribution (0) 
is no longer valid, i.e., the system enters a non-stationary 
nonequilibrium phase. 

The nature of the transition depends on the depth of 
the potential well. For go < §7 the potential well is deep 
enough to bind the interface to the wall at the transition 
point qc — 1, giving rise to a localized equilibrium distri- 
bution with a discontinuous transition. For q^ > 3, no 
localized solution exists at q — 1 and the transition be- 
comes continuous. The two transition lines are separated 
by a tricritical point at q^ = |, gc = 1- 

In order to prove the existence of the first-order line, we 
apply a transfer matrix formalism |3] . Defining a transfer 
matrix T acting in spatial direction by 



T, 
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go 
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(6) 



we compute the eigenvector (p oiT corresponding to the 
largest eigenvalue /x, which determines the steady-state 
properties of the system. For q = 1 the solution reads 



// = (z + l)/go, 
where h > I and 



Po = go, 



bh = z\ 



yi + 2qo - 3gg 1 
2(1 - go) 2 



(7) 



(8) 



The stationary density of exposed sites at the bottom 
layer is given by no = ^o/ShLo'^ft- ^^ ^^ nonzero for 
go < I and vanishes linearly at the tricritical point. This 
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FIG. 2. Phase diagram for p = 0.2. The inset shows the 
average time r for the interface to detach from the wall as 
a function of the system size A*' for p — qo = 0.2 inside and 
outside the coexistence regime. 



proves the existence of the first-order phase transition 
line in Fig. U 

In Ref. Ipl the density no and the interface width w — 



({h — {h))~y^^ at qo — q were found to scale as 



no ^ [qc - qf 



i<lc - qY 



(9) 



with the critical exponents xq = 1 and 7 = g. Using 
the transfer matrix approach, we can prove that these 
bulk exponents remain valid along the entire second or- 
der phase transition line, except for the tricritical point 
where xo = j = h- Moreover, approaching the tricritical 
point from the left along the first order transition line, it 
can be shown that the two quantities scale as 



«o~ iqo-<ioT"^ 

where Xq = 7* = 1. 



{Qo -Qo) ^ , 



(10) 



C. The case p ^ 1: In this case the dynamic rules 
do not satisfy detailed balance and the model cannot be 
solved using the previous methods. Performing Monte- 
Carlo simulations we determined the phase diagrams for 
various values of p. For p < 1 we find that the moving and 
the pinned phases coexist in a whole region of the parame- 
ter space rather than just on a line, as is the case for equi- 
librium first order transitions. As shown in Fig. H, the 
coexistence regime for p — 0.2 ends at the tricritical point 
q^ = 0.515(10), gc = 0.6868(2), where the second-order 
phase transition line starts. Unlike metastable states, the 
pinned phase is thermodynamically stable inside the co- 
existence regime, i.e., its life time r grows exponentially 
with the system size, as shown in the inset of Fig. 0. 

For p > 1, however, there is no region of phase co- 
existence and the phase diagram is similar to that of 




FIG. 3. Temporal evolution of a partially detached inter- 
face in the coexistence regime 5 = 0. 8, go =p = 0.2. Initially 
a large droplet is introduced by hand. The droplet quickly 
grows, reaches a triangular shape, and finally shrinks at con- 
stant velocity. 

Fig. |l|. For instance, for p = 2 the first-order phase 
transition line ends at the tricritical point q^ = 0.73(1), 
gc = 1.2326(3). For p ^ 1, we expect that Eqs. (9)-(10) 
still describe the scaling behavior at the tricritical point, 
although with different sets of critical exponents. 

In order to understand the mechanism leading to phase 
coexistence for p < 1, let us consider the evolution of a 
large droplet (an interval where the interface is detached 
from the bottom layer) in the vicinity of the upper termi- 
nal point of the coexistence regime g = 1, go = 0. Because 
of the RSOS constraint, the growing droplet eventually 
reaches an almost triangular shape with unit slope at 
the edges. The interface of the triangular droplet fluc- 
tuates predominantly by diffusion of pairs of sites with 
equal height. Inspecting the dynamic rules, it is easy 
to verify that these 'landings' of the staircase move up- 
wards with rate q and downwards with rate 1. Hence, 
for g > l,go = the droplet is stable with a life-time 
exponential in its lateral size. For g > 1 and go > 0, fluc- 
tuations of the bottom layer are biased to move upwards 
at the edges of the droplet. Thus the droplet grows and 
the interface eventually detaches from the bottom layer. 
On the other hand, if go = and gc < g < 1, fluctuations 
at the top of the triangular droplet are biased to diffuse 
downwards to the edges. Therefore, the droplet shrinks 
at constant velocity in a time proportional to its size, 
ensuring the stability of the pinned phase. 

As shown in Fig. 0, this robust mechanism for the elim- 
ination of droplets also works for positive values of go . If 
the interface detaches from the substrate over some dis- 
tance due to fluctuations, the resulting droplet grows and 
reaches an almost triangular shape. In the coexistence 
regime, the droplets are biased to shrink in a time propor- 
tional to its size, resulting in a stable pinned phase. How- 
ever, spontaneously created small islands next to the bot- 
tom layer contribute to the broadening of the droplets, 
reducing the range of phase coexistence. This explains 
why the upper boundary of the coexistence regime de- 
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FIG. 4. Schematic phase diagram in the p, g-plane, show- 
ing the regions of phase coexistence (PC) for fixed qo- Dis- 
continuous and continuous transition fines are represented by 
sofid and dashed fines, respectively. In the third panel the 
phase coexistence regime is enlarged artificially. The star de- 
notes the equilibrium transition point (see text). 

creases as go is increased. At the upper boundary the 
stationary density of exposed sites at the bottom layer 
no is found to change discontinuously. 

The transition hne above which the unbound phase is 
stable is independent of the growth rate qq. This line is 
the lower curve in Fig. |j, which is common to all four 
diagrams. For qq smaller than some threshold qo, the 
pinned and the unbound phases coexist in a certain re- 
gion of the phase diagram. As can be seen in the Figure, 
this region is bounded by two lines which intersect to the 
right at the equilibrium transition point p — q = 1. For 
Qa < qc,o — 0.399, this is the only intersection point of 
the two lines and the phase coexistence region extends 
down to p = 0. On the other hand, for go > qc,o the two 
lines also intersect on the left at another tricritical point, 
reducing the size of the region of phase coexistence. This 
region disappears at qo = Qo- On the basis of our numer- 
ical simulations, it is not possible to conclude whether go 
is equal or strictly smaller than 2/3. 

D. Discussion: Within a more general framework, 
the coexistence of the moving and the pinned phase may 
be viewed as follows. The evolution of the interface may 
be described in terms of the KPZ equation 



dth = DV^h + \{Vhf + C(a;, t) + V'{h) + 



«o 



(11) 



with positive heights h{x, t) > 0, where the velocity vo 
plays the role of g — gc . Clearly, for A = the transition 
takes place at vo = 0. For A > 0, the nonlinear term of 
Eq. (11) may be interpreted as an additional force acting 
on tilted parts of the interface in the direction of growth. 
This force supports the growth of droplets wherefore the 
interface detaches for any vo > 0. However, if A < 



this force acts against the direction of growth. Conse- 
quently, a sufficiently tilted interface does not propagate 
and may even move downwards. For wo > this leads 
to the formation of fluctuating droplets with a triangular 
shape and a finite slope at the edges. If the short-range 
force at the bottom layer is strong enough, such droplets, 
once formed, will shrink at constant velocity. Thus, the 
moving and the pinned phase can only coexist in those 
parts of the phase diagram where A is negative. In fact, 
as shown in , A is negative along the transition line for 
p < I and changes sign at p = 1. 

The phenomenon of phase coexistence was first ob- 
served in Toom's two-dimensional north-east-center vot- 
ing model 1^]. It was also shown that open boundaries 
in certain one-dimensional diffusive models may exhibit 
similar phenomena ||9[] . The model discussed in this work 
demonstrates that phase coexistence can also emerge in 
homogeneous one-dimensional driven systems. 
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